Shifts of optical frequency references based on spectral-hole burning in Eu 3+ :Y 2 Si05 
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Several properties of Eu 3+ :Y2SiOs spectral holes are measured, to assess the suitability of broad- 
band hole-patterns for use as laser- frequency references. We measure frequency shifts due to mag- 
netic fields, side-features of neighboring spectral holes, and changing optical probe power. A precise 
calibration of a temperature insensitive point is also performed, where the temperature-induced 
frequency shift is canceled to first order by the pressure-induced shift from the crystal's helium-gas 
environment. 



I. INTRODUCTION 

Frequency- stable laser-local-oscillators (LLO) are crit- 
ical components of the new generation of optical atomic 
clocks 1 ^—. The stability of such clocks is presently lim- 
ited by thermomechanical noise in the cavity-stabilized 
LLOs 4 - 5 . Quieter LLOs would enable optical clocks with 
improved stability, yielding faster comparison measure- 
ments among different clocks and more precise searches 
for the variation of fundamental constants 6-8 . Quieter 
LLOs may also yield lower-phase-noise microwave oscil- 
lators by frequency-division via femtosecond laser fre- 
quency combs^r— . 

Spectral-hole-burning (SHB) laser frequency 
stabilization^— in Eu 3+ :Y2SiOs has recently ap- 
proached the performance of state-of-the-art cavity- 
stabilized LLOsi 6 -, and indirect measurements suggest 
that the fundamental noise floor may be significantly 
lower. Furthermore, high-resolution spectroscopy 
has shown that the environmental sensitivity of 
Eu 3+ :Y 2 Si05 spectral holes to temperature, pressure, 
and accelerations are all smaller than that of Fabry- 
Perot cavities^ 6 -, and frequency drift rates can be lower 
than those of cavity-stabilized lasers 1 ^. These results 
suggest that spectral holes in Eu 3+ :Y 2 Si05 could form 
a frequency-stable flywheel for optical atomic clocks 
and potentially extend the coherence time of atom-laser 
interactions. 

Unlike Fabry-Perot cavities whose oscillation period 
is directly proportional to an artificial length, the opti- 
cal frequency of spectral holes depends primarily on the 
atomic Hamiltonian of dopant ions. Spectral holes are 
less sensitive to mechanical perturbations, which couple 
only weakly through the field of the host crystal, but 
more sensitive to ambient electric and magnetic fields. 
Moreover, probing of the spectral holes causes additional 
burning that can modify the hole shape and frequency. 
This work aims to determine the extent to which fluctu- 
ating atomic frequency shifts will limit the performance 
of SHB LLOs. 

In this work, we describe measurements of 
Eu 3+ :Y2Si05 spectral-hole properties including the 
magnetic field sensitivity, perturbations to the spectral- 
hole pattern due to side-holes and anti-holes, and 
frequency shifts due to fluctuations in the optical probe 
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FIG. 1. Experimental setup for measuring the properties 
of Eu 3+ :Y 2 Si0 5 spectral holes. A 580 nm dye laser is 
pre- stabilized to a spherical Fabry-Perot reference cavity. 
The laser is used to burn and probe spectral holes in two 
Eu 3+ :Y2Si05 crystals located in a single closed cycle cryostat. 
Spectroscopy is controlled by a field-programmable gate array 
(FPGA). AOM: acousto-optic modulator, EOM: electro-optic 
modulator, DET: detector. 



power. The current work is performed in a cryostat with 
independent control over temperature and pressure. 
Such environmental control enables a precise calibration 
of the temperature sensitivity, and it is shown that the 
linear dependence of the hole frequency on temperature 
can be eliminated. 



This paper proceeds as follows: Sec. |TT] de- 
scribes the experimental setup. Measurements of 
Eu 3+ :Y2Si05 spectral-hole properties are presented in 
Sees. IIIII (magnetic field sensitivity) and IIVI (tempera- 
ture and optical probe power sensitivity). Perturbations 
to spectral-hole patterns due to side-holes and anti-holes 
are discussed in Sec. [V] Finally, Sec. I VII summarizes the 
shifts of optical frequency references based on spectral- 
hole burning in Eu 3+ :Y2SiOs . 
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II. SETUP 



A schematic of the experimental setup is shown 
in Fig. Q] Spectral holes in cryogenically cooled 
Eu 3+ :Y 2 SiC>5 arc burned and probed by a dye laser 
that is pre-stabilized to a spherical Fabry-Perot reference 
cavit y 18 ' 19 . The laser frequency can be tuned coarsely 
by stepping among longitudinal modes of the reference 
cavity (3 GHz free spectral range) and finely by shift- 
ing the frequency of an acousto-optic modulator (AOM3) 
located between the laser and the cavity. The pre- 
stabilized laser has a fractional frequency stability of 
1.2 x 10 -15 between 0.5 s and 12 s and a typical drift rate 
of 10~ 14 /s. Eu 3+ :Y2Si05 exhibits two inhomogeneously 
broadened absorption features at 580.039 nm (site 1) and 
580.211 nm (site 2) that result from different crystal held 
shifts at the two distinct positions within the Y2SiC>5 unit 
cell where Eu 3+ can substitute for Y 3+ ; in this work we 
focus primarily on site 1 spectral holes because they dis- 
play a smaller temperature shift. 



All of the measurements presented in Sec. IIVI were 
performed on an identical pair of cylindrical samples 
of 1.0 atomic (at.) % Eu 3+ Y2Si05 with a diameter of 
10.0 mm and a length of 4.9 mm. For each sample the 
laser propagates parallel to the axis of the cylinder (b- 
axis of the crystal 2 ^), whose ends are wedged at a 2° 
angle to prevent etalon effects, and is polarized along the 
Dl crystal axis. The samples are cooled to 3.5 to 8.8 K 
in a closed-cycle cryostat that allows optical access to the 
cooled samples. Inside the cryostat, the two crystals are 
enclosed in a single hermetically sealed chamber that is 
filled with helium gas to provide independent control of 
temperature and pressure. The pressure is controlled by 
a thin tube that connects the sealed chamber to a room- 
temperature gas manifold and each crystal is mounted in 
a way that reduces the sensitivity to vibrations^. The 
measurements presented in Sees. Mil and [V] were made 
with samples of 0.1 and 0.5 at. % Eu 3+ :Y 2 Si05 in an 
optical flow cryostat described previously-^. All of the 
crystals used in this work were grown with a natural 
abundance of europium isotopes (47.8% 151 Eu and 52.2% 
153 Eu). 

Spectroscopy of Eu 3+ :Y2Si05 spectral holes is con- 
trolled by a field-programmable gate array (FPGA) that 
includes an embedded microprocessor. The laser beams 
(5 mm diameter) incident on each crystal can be tuned 
independently over a 600 MHz frequency range and 0.1 
to 20 /iW/cm 2 intensity by adjusting the radio-frequency 
(RF) drives of AOM 1 and 2. Phase modulation side- 
bands are imprinted on the lasers for PDH laser frequency 
stabilization^ by electro-optic modulators. Detectors II 
and 12 measure the power incident on the crystals, which 
is stabilized at the desired level by feedback to the am- 
plitudes of the RF drives of AOM 1 and AOM 2. Finally, 
detectors Tl and T2 measure the transmission of the 
laser through crystal 1 and crystal 2, respectively. 
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FIG. 2. The linear Zeeman effect in 

0.1 at. % Eu 3+ :Y 2 Si0 5 site 1 spectral holes at 4.5 K. 
(a) A magnetic field applied to spectral holes after they have 
been written splits the spectral hole into a symmetric set 
of spectral components. In this example a spectral hole is 
written in an applied magnetic field of G and measured 
in an applied magnetic field of 37 G along each of the 
crystal axes. The three measurements are offset vertically for 
clarity, (b) Frequency shift of the most prominent spectral 
components (indicated by vertical arrows in Fig. |2k) as a 
function of applied magnetic field along the crystal b axis. 



III. MAGNETIC FIELD EFFECTS 

The spectral-hole-burning mechanism in 
Eu 3+ :Y 2 Si05 is based on population redistribution 
among the hyperfine ground states. Therefore, changes 
in the magnetic field environment are expected to 
modify spectral-hole frequencies and lineshapes due to 
the Zeeman effect. However, the 7 Fn- 5 Dn transition 
of Eu 3+ has small magnetic moments in both the 
ground and excited states^, because only the nuclear 
angular-momentum is non-zero. 

Fig. [2^, shows the transmission spectrum of a hole writ- 
ten at zero applied field and subsequently measured in a 
magnetic field of 37 G (1 G = 0.1 mT) that is sequen- 
tially applied along each of the three crystal axes. Here 
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FIG. 3. The quadratic Zeeman shift in 

0.1 at. % Eu 3+ :Y 2 Si0 5 site 1 spectral holes at 4.5 K 
for magnetic fields applied along the b-axis, D2-axis, and 
Dl-axis. 



the spectral hole, which was a single peak at zero applied 
field, splits into a symmetric set of transmission peaks 
corresponding to transitions between pairs of Kramer's 
doublets. The Zeeman shifts are different for magnetic 
fields along the different crystal axes. Fig. [2}d shows the 
frequency shifts of the most prominent spectral compo- 
nents as a function of applied magnetic field along the 
crystal b axis. Because of the symmetry of the Kramer's 
doublets, the linear Zeeman shift averages to zero for 
each spectral hole, but it can contribute to broadening of 
the holes. The linear Zeeman shift of the most sensitive 
states is 3.8 kHz/G (magnetic field along the crystal b 
axis). To avoid broadening spectral holes at the 1 kHz 
FWHM level, the magnetic field fluctuations must be less 
than 0.1 G. 

Unlike the linear Zeeman effect, the quadratic Zeeman 
effect causes a frequency shift of spectral holes. Fig. [3] 
shows the area-weighted average (over the most promi- 
nent transmission peaks, indicated in Fig. [5^ by vertical 
arrows) frequency shift of a spectral hole written at zero 
applied field and measured in a nonzero applied mag- 
netic field. We measure a quadratic frequency sensitivity 
of 1.4(1.4) Hz/G 2 for magnetic fields along the crystal b 
axis, which is again the most sensitive direction. Note 
that we use the measured value as an estimate of the un- 
certainty because the atomic population of the original 
spectral hole is nearly twice the total population of the 
discernible Zeeman components. The remaining popu- 
lation is contained in spectral features whose center fre- 
quency is not readily discernible. The measured sensitiv- 
ity can be compared to the expected value of 0.64 Hz/G 2 
in free Eu 3+ atoms 23 ' 24 . In order to reach a fractional 
frequency instability of 10 -17 the magnetic field must be 
held at G with an instability better than 60 mG, or in 
a laboratory field of 0.5 G the magnetic field instability 
must be below 7 mG. 



IV. THERMAL EFFECTS 

At low temperature, the frequency shift (A/) of spec- 
tral holes in Eu 3+ :Y2Si05 exhibits a fourth-order tem- 
perature (T c ) dependence^. We measured Af = aT 4 
with a = 76(15) Hz/K 4 for site 1 and a = 250(50) Hz/K 4 
for site 2 over a temperature range of 2.5-5.5 K in a pre- 
viously described apparatus^. At 3.65 K, the slope is 
15 kHz/K for site 1 spectral holes and such a high sen- 
sitivity could make Eu 3+ :Y2Si05 crystals unsuitable for 
laser-frequency stabilization. 

Here we immerse the SHB crystal in a sealed helium- 
gas environment to mitigate the temperature shift 16 . 
Temperature increases cause the pressure of the helium 
gas to increase, which compresses the crystal. Because 
the spectral-hole frequency shift due to compression is 
negative and that due to a rise in temperature is posi- 
tive, the slopes of the two shifts can be made to cancel 
at a particular pressure of helium gas. Note, however, 
that this cancellation works only when the crystal is in 
thermal equilibrium with the surrounding gas. In our 
apparatus a tube connects the sealed crystal chamber to 
a room-temperature gas manifold through which helium 
can be added and removed. To suppress time-dependent 
frequency shifts associated with the nonzero thermaliza- 
tion time of the two-reservoir system, we have minimized 
the volume of the room-temperature portion such that 
it is approximately equal to the volume of the cryogenic 
portion. In the limit V r /V c -C T r /T c where V r (V c ) and 
T r (T c ) are the volume and temperature of the room- 
temperature (cryogenic) portion of the helium gas, the 
pressure at which first order temperature shifts cancel is 
Pa = —AaTy/P, where To is the operating temperature 
and the spectral-hole pressure shift has been measured 16 
to be P = -211.4 Hz/Pa for site 1 and /3 = -52.0 Hz/Pa 
for site 2. The residual quadratic temperature sensitivity 
is A/ = 6aT 2 (T c - T ) 2 

Figures 3k and HJj show the temperature-dependent 
frequency shift in 270 Pa of helium gas. The frequency 
difference between the pre-stabilization Fabry-Perot cav- 
ity and a laser that is locked to a pattern of spectral 
holes is recorded as a function of time while the tem- 
perature is changed in discrete steps near the insensitive 
point. Figure |4p shows the frequency shift as a func- 
tion of the crystal temperature. The residual quadratic 
temperature shift at 3.65 K for site 1 spectral holes is 
6.5 kHz/K 2 , which implies that in order to reach 10~ 17 
fractional frequency instability, the temperature instabil- 
ity needs to be below 1 mK. If the temperature is 5 mK 
away from the insensitive temperature, a stability better 
than 80 fiK is required. 

Figure 0p shows the temperature stability of the crys- 
tal chamber inside the closed-cycle cryostat whose cold 
head is thermally stabilized by a feedback loop. The sta- 
bility shown here is measured by an out-of-loop sensor 
located on the copper crystal chamber. The chamber is 
thermally decoupled from the cold head by brass washers, 
which improves the temperature stability of the chamber 
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FIG. 4. Calibration of the temperature insensitive point of 
1.0 at. % Eu 3+ :Y 2 Si0 5 site 1 spectral holes at 270 Pa. (a) 
Time-trace of the relative frequency between spectral-hole 
patterns in crystals 1 and 2 and the pre-stabilization refer- 
ence cavity. The temperature is stepped in the vicinity of 
the insensitive temperature causing shifts in the frequency of 
the hole pattern, (b) Fits of the frequency shifts reveal both 
the insensitive temperature and the quadratic temperature 
sensitivity, (c) Temperature stability (Allan deviation) of the 
crystal chamber. The total measurement duration is 2 x 10 s. 



at short times but degrades it at long times. Tempera- 
ture stability is better than 1 mK for the entire range of 
measurement times and better than 80 /iK for durations 
between 250 s and 50000 s. 

Local heating of the crystal by the laser may cause an 
additional thermal frequency shift. However, this effect 
is reduced because the helium gas inside the chamber 
thermalizes the crystal with the chamber walls. We have 
measured an upper bound for the spectral-hole frequency 
shift due to a step change in laser power. We locked two 
beamlines of the laser separately to a pattern of site 1 
spectral holes in each crystal at 0.4 /iW optical power for 
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FIG. 5. Side-hole and anti-hole spectrum due to the 
Eu 3+ :Y2SiOs site 1 spectral-hole hyperfine structure as a 
function of spectral-hole width and Eu 3+ doping. In each 
transmission measurement, a single spectral hole is burned to 
a peak amplitude of roughly 600 in the arbitrary transmis- 
sion units of the figure. The legend labels the transmission 
measurements by specifying the Eu 3+ doping level, as well as 
the crystal temperature and the FWHM of the spectral hole. 
The amplitude of the side-holes and anti-holes is reduced by 
writing narrower spectral holes because the side-holes and 
anti-holes exhibit inhomogeneous broadening. 

400 s, then increased the laser power incident on crystal 
2 to 1.3 jLtW for 50 s (leaving the laser power incident 
on crystal 1 unchanged). The relative frequency shift 
was less than 1 Hz, corresponding to an optical power 
sensitivity of less than 1 Hz//iW. Current power stability 
is 4 nW, so laser power fluctuations are unlikely to limit 
the performance of this system at the 10 -17 fractional 
frequency stability level. 

V. SIDE-HOLES AND ANTI-HOLES 

Side-holes and anti-holes™ perturb broadband 
spectral-hole patterns, which may be the best way to 
implement a SHB LLO that can operate continuously. 
Figure [5] shows the side- hole and anti-hole spectrum of a 
single spectral hole for several spectral-hole widths and 
two Eu 3+ doping levels. Inhomogeneous broadening in 
the Eu 3+ :Y2Si05 hyperfine structure leads to side-holes 
and anti-holes that are broader than the main spectral 
hole feature. As a result, the amplitude of the side-holes 
and anti-holes can be reduced by writing narrow spectral 
holes. As the Eu 3+ doping level increases, the inho- 
mogeneous broadening in the Eu 3+ :Y2Si05 hyperfine 
structure increases. Figure [5] shows the widths and 
amplitudes of the side-holes and anti-holes of a single 
spectral hole for two Eu 3+ doping levels. Increasing the 
doping from 0.1 % to 0.5 % results in a mean increase 
in the widths of the side-holes and anti-holes by a factor 
of 3.8. Perturbations from side-holes and anti-holes to 
broadband patterns can therefore be reduced by choos- 
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FIG. 6. FWHM and amplitude of the side-holes (positive) and 
anti-holes (negative) due to the Eu 3+ :Y2SiOs site 1 spectral- 
hole hyperfine structure as a function of Eu 3+ doping. In each 
transmission measurement, a single 11.6 kHz (0.1% doping at 
7.3 K) or 83 kHz (0.5% doping at 8.8 K) FWHM spectral hole 
was burned to full transparency; the transmission amplitude 
is normalized to the transmission of the main hole. 



TABLE I. Properties and systematic frequency shifts of the 
Eu 3+ :Y 2 Si0 5 crystal. 



Parameter 



Site 1 



Site 2 



Units 



Wavelength 
Bulk modulus 2 ^ 
Volume^ 
Magnetic field 
Electric neld 2 ^ 
Temperature^ 
Pressure^ 
Optical power 
Acceleration^ 



580.039 

134.8 
0.055 
1.4(1.4) 

< 0.2 
6.2 

-211.4(4) 

< 1 

7 x 10~ 12 



580.211 
134.8 
0.014 



21 

-52.0(7) 



GPa 

(Af/f)/(AV/V) 

Hz/G 2 

Hz/(V/m) 2 

kHz/K 2 

Hz/Pa 

Hz//jW 

(A///)/ff 



a at 3.7 K and 270 Pa (3600 Pa) helium gas pressure for site 1 (2) 

ing the highest dopant concentration that is compatible 
with the requirements for spectral-hole lifetime and 
coherence. 

The effect of side-holes and anti-holes on the perfor- 
mance of a SHB LLO can be estimated by considering 
how much the center frequency of a spectral hole will be 
shifted by an overlapping side-hole or anti-hole. If two 
2 kHz FWHM spectral holes of equal depth are burned in 
0.5 at. % Eu 3+ :Y 2 Si05 with worst-case frequency spac- 
ing, then side-features of the second hole will shift the 
center frequency of the first hole by only 2 mHz. Thus, 
SHB LLO frequency noise due to side-holes and anti- 
holes can be less than 10 -17 . 



tions in magnetic field, temperature, pressure, and opti- 
cal power to be below 10~ 17 in our experimental setup. 

While the linear Stark coefficient of Eu 3+ :Y 2 Si05 site 
1 spectral holes has been measured^ to be 
350(50) Hz/(V/m), the more relevant quadratic 
Stark effect has not yet been measured. Similar to the 
Zeeman effect, when the amplitude of the electric field 
is small, the linear Stark effect causes spectral-hole 
broadening and the quadratic Stark effect causes the 
center frequency of the spectral holes to shift. The 
linear Stark effect data of Graf— are consistent with a 
zero quadratic Stark effect with a large uncertainty of 
±0.2 Hz/(V/m) 2 for Eu 3+ :Y 2 Si0 5 site 1 spectral holes. 

The rate of spectral-hole frequency drift due to back- 
ground electric and magnetic field drift can be estimated 
using the linear and quadratic field sensitivities together 
with a measurement of the rate of spectral-hole broad- 
ening. We have measured the rate of Eu 3+ :Y 2 SiOs site 
1 spectral-hole broadening in the absence of broadening 
due to continuous probing by burning a 2.1 kHz FWHM 
spectral hole and measuring the width every 30 minutes. 
The observed growth rate of the spectral-hole FWHM, 
0.04 Hz/s, indicates that the rate of change of background 
fields is likely below 6 x 10~ 5 (V/m)/s and 5 x 10" 6 G/s. 
In a laboratory magnetic field of 0.5 G, such a magnetic 
field drift would cause less than 7 x 10~ 6 Hz/s of spectral- 
hole frequency drift. 



VII. CONCLUSION 

Spectral-hole burning in Eu 3+ :Y2Si05 shows much 
promise for improved LLOs. The systematic frequency 
shifts of Eu 3+ :Y 2 Si05 spectral holes due to magnetic 
field, temperature, helium gas pressure, and optical probe 
power fluctuations are all small enough to allow laser 
frequency stabilization at the 10 -17 level in the current 
apparatus. 

We have shown qualitatively that perturbations to the 
spectral-hole pattern due to side-holes and anti-holes 
can be suppressed by using narrower spectral holes and 
higher doping levels, but more work is required to deter- 
mine the contribution of side-holes and anti-holes to SHB 
LLO frequency noise. Further studies of the quadratic 
Stark shift and the thermomechanical noise floor are also 
needed, but we are optimistic that spectral-hole burning 
can reduce the frequency noise of lasers and extend their 
coherence time. 
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